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ABSTRACT 

We report the kinematic discovery of the Pisces Stellar Stream (PSS), at Galactic 
longitude I « 135° and —39° < b < —36°. We originally identified this halo substructure 
from velocities of red giant branch stars in Sloan Digital Sky Survey (SDSS) Data 
Release 8, and confirmed its presence in turnoff stars from SDSS photometric data. The 
PSS is a narrow, kinematically cold tidal stream, with cr v ,o ~ 8 km s _1 . Its metallicity is 
[Fe/H] ~ —2.2, with ~ 0.3-dex dispersion. The color-magnitude signature of the stream 
turnoff, combined with our measured metallicity, places the PSS at a distance of 35±3 
kpc. The Pisces Stellar Stream is the same as the previously announced "Triangulum 
stream" and part of the proposed "stream a". We rule out an association of the PSS 
with other previously known Milky Way substructures in the same region of the sky. 

Subject headings: Galaxy: structure - Galaxy: kinematics and dynamics - Galaxy: 
stellar content 



1. Introduction 



Many stellar streams and substructures resulting from the tidal disruption of satellite dwarf 
galaxies and star clusters under the gravitational influenc e of the Milky Way have been identified 



and mapped using deep, large-area photometric surveys 
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2). Streams and cloud- like substructures have also been ob- 



served kinematically via spectroscopic surveys of tracers such as RR Lyrae, M-giants, or turnoff 
stars (e.g., Vivas et al. 2001; Duffau et al. 2006; Carlin et al. 2012; Majewski et al. 2012; Sheffield 
et al. 2012). 



Here we report a stellar tidal stream discovered kinematically, which (based on its location in 
the constellation Pisces) we dub the "Pisces Stellar Stream" (PSS). This stream was identified as 
a grouping of red giant branch stars at a common distance and radial velocity during the course 
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Fig. 1. — Line of sight veloc ity in the Galactic standard of rest frame, v gsr (calculated using 



Equation 1 from lLi et al.ll2012l ). as a function of Galactic latitude, for stars from SDSS DR8. These 
are selected as metal-poor giant stars in the south Galactic cap, with 120° < I < 160°. Two 
separate bands of stars with v gsr < km s _1 sweep across the lower portion of this figure, moving 
to lower v gsr with increasing latitude. These two structures are associated with the Cetus Polar 
Stream and the Sagittarius tidal stream. An additional clump of stars is seen at b ~ —37° and 
fg SI . ~ 90 km s ; we show in this work that these stars are part of a stellar stream at a distance 



of ~ 35 kpc. 



of a spectroscopic study of the Cetus Polar Stream (CPS). In the CPS study (Yam et al. 2013), 
we selected metal-poor (—3.0 <[Fe/H]< —1.5) stars with surface gravities consistent with being 
giants (0 < log g < 4.0) from the SDSS DR8 spectroscopic data in the south Galactic hemisphere. 
The Cetus stream was clearly present in velocity versus Galactic latitude plots of these stars, along 
with an additional feature consisting of a few stars clumped at Vg Sr ^ 90 km s ^ . Upon further 
exploration, these stars turn out to nearly all be drawn from the same small region of sky and 
have similar magnitudes (and thus similar distances). In this Letter, we show evidence that this 
grouping of stars is part of a stellar tidal stream at a distance of roughly 35 kpc from the Sun in 
the constellation Pisces. This stream is a piece of "stream a" in the matched-filtered density maps 
of Grillmair (2012), and is the same as the "Triangulum stream" in Bonaca et al. (2012). We chose 
to rename the stream because our detection is within the Pisces constellation, the stream does not 
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pass through Triangulum, and we cannot verify whether the stream extends the length of "stream 
a." This Letter presents the first measurement of velocities and metallicities of stars in this stream. 

2. Observation and Data Analysis 

In Figured] we plot v gST vs. Galactic latitude for stars between Galactic longitudes of 120° — 
160°, < logg < 4 (giant stars), -2.8 < [Fe/H] < -1.8, 0.4 < (g - r) < 0.9 and g > 15.5. 
The color and magnitudes were chosen to exclude both turnoff stars whose surface gravity was 
miscalculated and evolved post-RGB stars, while retaining likely red giant branch stars. The most 
prominent features in Figure Q] are swaths of RGB stars at v gsr < km s _1 , running across most 
of this plot, which are associated with the Sagittarius and Cetus streams (Newberg et al. 2009, 
Koposov et al. 2012, Yam et al. 2013). An additional overdensity of metal-poor giant stars can be 
seen around b = —37° and v gST = 90 km s . 

We now show that this velocity clump represents a real structure in the Galactic halo. From 
the sample of metal-poor red giant branch (RGB) stars shown in Figure [H we selected a smaller 
velocity range of 75 km s^ 1 < v gST < 115 km s^ 1 . We plotted the positions of stars in Galactic 
coordinates, and noted that many of these velocity-selected giant stars with metallicities between 
—2.8 <[Fe/H]< —1.8 are concentrated in a small fraction of a single SDSS spectroscopic plate 
centered at I = 136° and b = —35° (see Figure [2]). 

Figured was produced using the matched- filter technique described in detail by Grillmair (2009), 
which is in turn a refinement of the Rockosi et al. (2002) method. This method can be used to 
create weighted density maps of stars from within a given color-magnitude filter (based on either 
an empirical ridgeline or a theoretical isochrone) placed at a given distance; the weights empha- 
size regions of the filter expected to contain the highest ratio of stream stars to foreground stars. 
The map in Figure [2] shows filtered star counts from SDSS DR8 using a filter based on a Padova 
isochrone (Marigo et al. 2008, Girardi et al. 2010) with Z = 0.0001 ([Fe/H]~ -2.28), age 12.7 Gyr, 
at a distance of 33.6 kpc (note that the filter's response will be highest near the main sequence 
turnoff of the filter, at colors where foreground contamination is minimal). The most prominent 
feature (the blob at a, 5 ~ 23.5°, 30.5°) is due to the M33 galaxy. A narrow stream is visible 
running nearly vertically across this plot, from (a, 5) ~ (24°, 23°) to (a, 5) ~ (21°, 35°). 

We selected a region of sky based on the positions of these velocity-selected, metal-poor giants 
that we call the "on-stream field." This is defined as the area within ±0.75° in right ascension 
about the line 5 = — 4.33a + 126.46°, confined to 22° < 5 < 28°; the linear trend was fit to apparent 
density peaks in the matched- filter star count map (Figure [2]). The width of the on-stream field 
was chosen to include all of the spectroscopic candidates we originally identified in this region, and 
its length covers the portion of the sky where the stream is obvious in the filtered star counts map. 

An off-stream region was defined that contains an area eight times larger than that of the on- 
stream field. This comparison region was made up of two areas (seen as dashed boxes on Figure [2]) 
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Fig. 2. — Left panel: map of filtered star counts in the region surrounding the Pisces Stellar 
Stream; darker shades represent higher densities. This map was created by extracting stars from 
SDSS DR8 within an isochrone filter of age 12.7 Gyr, metallicity Z = 0.0001 ([Fe/H]« -2.28), 
and placed at a distance of 33.6 kpc, using the matched filter technique of Grillmair (2009) with 
0.1° bins. The image was smoothed with a Gaussian kernel three pixels wide. The large "blob" 
at (a, 5) ~ (23.5°, 30.5°) is M33. Passing near M33 on the sky is the PSS, which is visible as 
a narrow feature running diagonally from (a, 8) ~ (24°, 23°) to (21°, 35°), and perhaps further. 
Red lines and dashed lines show the 'on-stream' and 'off-stream' regions of the sky, respectively. 
Note that the broad feature at a > 24° and 19 < 5 < 22° is a data artifact that arises because 
of poor seeing when this particular SDSS stripe was observed. Right panel: Positions of SDSS 
spectra of metal-poor (—2.8 <[Fe/H] < —1.8) giant stars (small open triangles). Filled red dots 
indicate stars with velocities between 75 km s _1 < v gsr < 115 km s _1 . The SDSS plate centered 
at (a, <5) ~ (24.7°, 23.7°) contains most of the velocity-selected red points. Within the plate, these 
metal-poor, velocity-selected giant stars are mostly clustered in a narrow line of 11 stars. 



defined by: 



- 0.231(5 + 21.455° < a < -0.231(5 + 27.455°, 22° < 6 < 28°; 



(1) 
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- 0.2315 + 30.956° < a < -0.2315 + 36.956°, 22° < 5 < 28°. (2) 

In the top left and top right panels of Figure [3] we show color magnitude diagrams (CMDs) of 
all stars having SDSS DR8 spectra within the on- and off-stream fields, respectively. We also show 
subsets of the spectra with properties expected for the PSS. The subscripts on the magnitudes go 
and {g — r)o mean that these magnitudes have been extinction corrected based on the Schlegel et 
al. (1998) maps (using the extinction values provided for each star in the SDSS database). 

The metal-poor giant stars selected using the velocity criteria (i.e., red points in the top 
panels of Figure [3]) are mostly clustered around [(g — r)o,<?o] = [0.55, 18.8] in the on-stream field 
(top left panel of Figure [3]), while the more metal-rich stars (light blue points) scatter throughout 
the diagram. We suggest that these stars are part of the red giant branch of a single structure. 
The comparison off-stream field (top right panel) contains several stars that match our velocity, 
metallicity, and giant-selection criteria, but don't have a coherent distribution that could reasonably 
correspond to a giant branch. In the on-stream field there are a total of 448 stars with spectra 
in this diagram: 32 light blue points (velocity-selected giants) and 11 red points (velocity and 
metallicity-selected giants). The off-stream area has 1584 total spectra, with 90 light blue points 
and 14 red points. To calculate potential foreground contamination we selected stars with go < 19.2 
that are within 0.1 magnitudes in go, (d ~ r )o of the Padova isochrone with age 12.7 Gyr, [Fe/H] 
= -2.3, at a distance of 33 kpc. There are 403 and 124 spectra within this filter in the off- and 
on-stream regions, respectively; a factor of 3.25 times as many spectra were observed in the off- 
stream region. Among the off-stream spectra in this isochrone filter, there are only 2 metal-poor 
(-2.8 < [Fe/H] < -1.8) red giants (log g < 4.0) with 75 < v gsr < 115 km s _1 . Scaling by the factor 
of 3.25, this leads to an expectation of ~ 1 background "interloper" in our on-stream sample. Our 
finding of 10 candidate stream stars in this region (satisfying the same criteria) is thus a significant 
excess relative to neighboring regions. Furthermore, these stars are clustered in a small region of 
the CMD where we expect RGB stars. One additional candidate is a blue star that is consistent 
with being a horizontal branch star at the same distance. Table 1 shows the properties of the eleven 
stars classified as PSS candidates. 

On the top left panel of Figure [3] we overplot three isochrones shifted to distances (28, 33, 
and 38 kpc) that roughly match the RGB, with a metallicity of Z = 0.0001 ([Fe/H]~ -2.28) and 
age 12.7 Gyr. Our choices of isochrone parameters will become clear later in this paper, as these 
were motivated by the main sequence turnoff (MSTO) visible in our background-subtracted Hess 
diagram, our spectroscopic metallicity measurement, and examination of our matched-filter star 
count map. For now, we note that the RGBs of the fainter two isochrones match the clump of red 
points in Figure [3] well. 

The surface gravity cut selects asymptotic giant branch (AGB) and blue horizontal branch 
(BHB) stars in addition to the RGB stars we are interested in finding. To further isolate RGB 
stars, we chose only stars with color (g — r)o > 0.4 from the velocity-selected sample shown in 
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Fig. 3. — Top row: color magnitude diagrams of stars within the on-stream (upper left panel) and off-stream (upper right 
panel) fields defined in Section 2. Small black dots represent all stars with SDSS DR8 spectra. Light blue points in both 
panels represent giant stars with Galactoccntric line-of-sight velocities 75 km s — 1 < v gBI < fl5 km s . Red points have 
an additional metallicity constraint of —2.8 < [Fc/H] < —1.8. The left (on-stream) panel has a clear concentration of stars 
grouped at 0.5 < (g — r)o < 0.7, go ~ 18.8, indicative of a red giant branch. Overplotted Padova isochrones correspond to a 
metallicity of [Fe/H] —2.3 and an age of 12.7 Gyr. Blue, green, and red isochrones correspond to distances of 28, 33, and 38 
kpc, respectively. Bottom row: Color magnitude Hess diagrams showing stellar density in the on-stream field after subtraction 
of the local background (off-stream, scaled down by a factor of eight) contribution. The left panel shows the residuals; the 
right panel is the same Hess diagram with Padova isochrones from the upper left panel overplotted. Darker pixels signify excess 
on-stream stars. A main sequence turnoff is apparent at go ~ 21 in the lower panels. 

Figure [3] (light blue points), and histogrammed their metallicities in the upper panel of Figure [H 
Because the off-stream field has many more stars than the on-stream selection, we normalized 
the off-stream counts by the relative number of stars with [Fe/H]> —1.5 (43 in the off-stream vs. 
11 in the on-stream region). The distributions of velocity-selected RGB candidate stars in the 
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on- and off-stream fields are clearly different. The on-stream field has a metal-rich component at 
—1.0 < [Fe/H] < —0.4 that is consistent with the off-stream (i.e., background) stars, but also a 
different, more metal poor, population at [Fe/H] < —1.8. Also shown in the left panel of Figured] 
is a Gaussian that we fit to the residuals from subtracting the scaled off-stream counts from the 
on-stream data; this yields ([Fe/H]) = —2.24, with <7r Fe / H ] = 0.32. Because this low-metallicity 
peak is not present in the off-stream fields, we identify this as the PSS metallicity signature; this 
is the reason for the metallicity selection of —2.8 < [Fe/H] < —1.8 for the red points in Figure EJ 

We performed a similar procedure to find the velocity of stream candidates. The right panel 
of Figure H] shows velocities of metal-poor ([Fe/H]< —1.5) red giants, with the off-stream number 
counts scaled down by a factor of 3.16 so that the number of red giants with [Fe/H]> —1.5 is equal 
in the two samples. Again, we subtracted the scaled off-stream histogram from the on-stream 
data, and fit a Gaussian to the residuals. The fit (shown as a dotted curve in Figure U|) gives 
(v gsr ) = 90.7 km s" 1 , with a v = 8.7 km s _1 . 

Having confirmed the presence of a population of stars having common velocity and metallicity, 
we next examine the photometric data in this region of the sky to look for additional evidence of 
the substructure. We selected all stars with photometric data from the SDSS DR8 database, within 
the on- and off-stream fields described above. We created CMD Hess diagrams for both the on- and 
off-stream fields using extinction-corrected go magnitudes and (g — r)o colors. The Hess diagram 
for the off-stream field was normalized by dividing number counts in each bin by 8 to account 
for the 8x larger area relative to the on-stream region. We subtracted the normalized off-stream 
Hess diagram from that of the on-stream field to create a difference plot. The residuals from this 
subtraction can be seen in the bottom left panel of Figure [3l 

Figure [3] (lower panels) reveals an overdensity of stream stars between 0.2 < (g — r)o < 0.4 and 
20.5 < go < 21.6, consistent with a population of F-turnoff stars from the stream whose velocity 
signature we have found in RGB stars. The turnoff magnitude is clearly fainter than that of the 28 
kpc isochrone, and appears more likely to be between the 33-kpc and 38-kpc isochrones. Since the 
stars are in a narrow line across the sky, we surmise that the turnoff stars and spectroscopically 
selected RGB stars are part of an old, metal-poor tidal stream at a distance of ~ 35 ± 3 kpc. 

3. Discussion and Conclusion 

We have identified the kinematical and chemical signature of a likely stellar tidal stream 
based on spectra in an SDSS plate located at / rz> 135° and between —39° < b < —36°. Stars in 
this structure, which we dub the Pisces Stellar Stream (PSS), have line-of-sight velocities in the 
Galactic standard of rest frame between 77 km s _1 < i> gsr < 114 km s" 1 . We estimate the mean 
velocity for the 11 stars using a maximum likelihood method (e.g., Pryor & Meylan 1993; Kleyna 
et al. 2002), and find (v gsI ) = 95.6 ±3.1 km s _1 and an intrinsic velocity dispersion (accounting for 
individual velocity errors; the typical error for an individual star is ~ 7.0 km s _1 ) of <7 v> o = 7.6±3.0 
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Fig. 4. — Left panel: Histogram showing the samples that were depicted as colored (blue and 
red) points in Figure El with an additional cut requiring stars to have (g — r)o > 0.4. Grey filled 
bars represent the on-stream sample and grey hashed bars represent stars in the off-stream field. 
The off-stream bins have been normalized so that the total number of stars with [Fe/H]> —1.5 is 
equal in the on- and off-stream samples. The on-stream stars consist of two different populations: 
the metal-rich peak centered on [Fe/H] ~ —0.7 roughly matches the "background" stars of the 
off-stream field, while a more metal-poor population is evident between —2.8 < [Fe/H] < —1.8 
only in the on-stream data. We fit a Gaussian to the residuals from subtracting the off-stream 
("background") from the on-stream histogram; the result is shown as a dotted curve centered at 
[Fe/H]=-2.24, with o\F e /m = 0.32. The sum of the off-stream counts and the Gaussian fit (dashed 
black line) matches the on-stream data well. Right panel: GSR-frame velocities of metal-poor 
([Fe/H]< —1.5) red giant stars in the on- and off-stream fields. Line styles, fill styles and scaling 
are as in the left panel. As in the metallicity plot, we subtracted the off-stream from the on-stream 
and fit a Gaussian to the residuals, this yields (v gsr ) = 90.7 km s _1 , with a v = 8.7 km s . The 
apparent peak located at v gsr ~ 50 km s -1 of stars appears to be a different population from the 
PSS. Three of the stars are separated by more than a magnitude from the PSS locus on the CMD 
and the final star has a velocity more than 6a from the measured value of the stream. 



km s . The metallicity of these 11 stars, determined from a Gaussian fit to the residuals after 
subtracting an off-stream field, is ([Fe/H]) = —2.24 with a dispersion of <7[F pm] = 0-32, wh ich is 



only slightly larger than the SDSS metallicity measurement error of ~ 0.2 dex (jLee et al.ll2008l ). We 
isolate the PSS in both a background-subtracted color-magnitude Hess diagram and a matched- 
filter star count map. From these, we derive a distance to the stream of 35 ± 3 kpc. No distance 
gradient was apparent over the short stretch of the stream we isolated. The mean position of the 
11 spectroscopically confirmed stars is (a, 5) = (23.82°, 23.53°). 

The Pisces stream is a small segment of the structure identified as "stream a" by Grill- 
mair (2012) in a matched-filter star count map of the south Galactic cap. This structure has 
also been reported by Bonaca et al. (2012) based on star count maps using a similar technique. 
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Bonaca et al. found a distance of 26 ± 4 kpc to the stream, which they dubbed the "Triangu- 
lum stream". However, this distance was based on a metallicity of [Fe/H] = -1.0 estimated from 
isochrones; our measured spectroscopic metallicity of [Fe/H] = -2.2 requires a more distant, ~ 35 
kpc isochrone to match the stream in the CMD. Our finding that the stream is kinematically cold, 
narrow in width, and essentially mono-metallic, confirms the suggestions of both Grillmair (2012) 
and Bonaca et al. (2012) that this stream is the remnant of a globular cluster, or possibly an 
ultrafaint dwarf galaxy. 

This stream is distinct in distance and line-of-sight velocity from other known substructures 
in this region of the sky. The Triangulum- Andromeda (TriAnd) stellar structure (Rocha-Pinto et 
al. 2004) is a large, diffuse substructure. TriAnd is ~16-25 kpc from the Sun (Majewski et al. 2004), 
much closer than our 35 kpc estimate for the PSS. Moreover, our mean velocity of (v gsr ) = 95.6 
km s -1 at I ~ 135° is much higher than the ~ 30 km s _1 found by Rocha-Pinto et al. (2004) for 
TriAnd at this position. We thus consider an association between the Pisces stream and TriAnd 
to be unlikely. The PSS is also unlikely to be associated with the CPS, as the CPS has v gST ~ —60 
km s _1 at b ~ —38°, compared to our ~ 90 km s _1 velocity for the Pisces stream at this position. 

We thank the anonymous referee for insightful comments. This work was supported by NSF 
grants AST 09-37523, AST 10-09670, and AST 12-39904. Funding for SDSS-III has been pro- 
vided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science 
Foundation, and the U.S. Department of Energy Office of Science. The SDSS-III web site is 
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Table 1. Pisces Stellar Stream Candidates 



ra (deg) dec (deg) V gsr (km s x ) [Fc/H] log .9 # (0 - r)o (.9 - i)o (u ~ 9)0 



24.18291 


22.93643 


113.5 


-2.54 


2.30 


18.14 


0.59 


0.84 


1.23 


24.21572 


22.95978 


93.1 


-1.87 


2.92 


18.98 


0.53 


0.73 


1.15 


23.80551 


23.47830 


96.7 


-2.37 


2.72 


18.89 


0.50 


0.67 


0.98 


23.92007 


23.39030 


85.5 


-2.13 


2.18 


18.71 


0.52 


0.73 


1.22 


23.24328 


23.19339 


80.9 


-2.12 


3.70 


17.84 


-0.14 


-0.26 


1.28 


23.57924 


24.39110 


101.7 


-2.65 


1.69 


19.12 


0.56 


0.76 


1.21 


23.41445 


24.34511 


93.1 


-2.07 


3.09 


19.10 


0.52 


0.72 


1.12 


23.78164 


23.88001 


77.4 


-2.38 


1.85 
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